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ABSTRACT 


A  simple  in— phase/quadrature  (I/Q)  demodulator  architecture  which  can 
measure  the  amplitude,  phase  and  instantaneous  frequency  of  radar  signals  is  examined  in 
this  report.  This  architecture  has  a  potential  for  radar  ESM  applications  where  wide 
instantaneous  bandwidth  and  simple  algorithms  for  extracting  the  modulation 
characteristics  of  radar  signals  are  required.  This  I/Q  architecture  can  meet  the 
requirement  by  splitting  an  incoming  signal  into  its  in— phase  and  quadrature  components 
using  analog  hardware.  However  in  practice,  there  are  amplitude  and  phase  imbalances 
between  the  two  components  and  DC  offsets,  which  can  introduce  large  systematic  errors  in 
the  measurements.  In  this  report,  we  present  novel  techniques  which  can  greatly  reduce  the 
systematic  errors  and  improve  the  accuracy  of  the  measurement.  A  time-domain  analysis 
of  the  systematic  errors  is  given.  A  calibration  technique  which  can  be  used  to  correct  for 
the  imbalances  and  o&ets  is  discussed.  The  effect  of  noise  on  the  accuracy  of  the 
measurement  is  also  examined.  Imbalance  errors  and  DC  offsets  of  I/Q  networks  are 
measured  and  analyzed.  Finally,  a  post— processing  technique  emplo^ng  moving  averages, 
which  is  shown  to  be  effective  for  improving  the  output  signal— to-noise  ratio  and  reducing 
systematic  errors,  is  also  presented. 


RESUME 

Une  architecture  de  d4modulateur  i  quadrature  de  phase  pour  la  mesure  de 
I'amplitude,  de  la  phase  et  de  la  frequence  instantande  des  si^aux  radar  est  exazninee  dans 
ce  rapport.  Cette  architecture  est  prometteuse  pour  les  applications  de  "mesures  de 
soutien  dlectronique  radar"  ou  une  grande  largeur  de  bande  instantanee  et  des  algorithmes 
simples  permettant  d’extraire  la  modulation  des  signaux  radar  sont  requis.  Cette 
architecture  rencontre  ces  exigences  en  sdndant,  avec  des  composants  analogiques,  un 
signal  d’entrde  en  ses  deux  composantes  en  quadrature  de  phase.  Cependant  des 
debalancements  d’amplitude  et  de  phase  entre  les  deux  composantes,  ainsi  que  des  biais  de 
tension,  peuvent  en  pratique  introduire  des  erreurs  systematiques  assez  grandes  dans  les 
mesures.  Une  analyse  temporelle  de  ces  erreurs  systematiques  est  presentee.  Une 
technique  de  calibration  permettant  de  corriger  les  debalancements  et  les  biais  de  tension 
est  discutee.  L’effet  du  bruit  sur  la  precision  des  mesures  est  aussi  examine.  Les  erreurs 
dues  aux  debalancements  et  aux  biais  de  tension  dans  les  demodulateurs  a  quadrature  sont 
mesurdes  et  analysees.  Enfin,  une  technique  de  filtrage  qui  permet  de  rdduire  les  erreurs 
systematiques  est  presentee. 
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EXECUTIVE  SUMMARY 


Due  to  the  increasing  density  and  complexity  of  radar  signal  waveforms,  it  is  no 
longer  easily  possible  to  sort  out  uniquely  and  identify  each  radar  emitter  using 
convention^  signal  parameters  such  as  pulse  width,  radio  frequency  (RF),  amplitude  and 
pulse  repetition  frequency.  As  a  result  there  is  a  requirement  for  a  radar  Electronic  Support 
Measures  (ESM)  receiver  to  measure  the  modulation  characteristics  of  radar  signals  and  to 
provide  additional  parameters  which  can  be  used  to  identify  unambiguously  eaw  type  of 
radar  emitter.  With  the  advent  of  fast  A/D  converters  and  high-sp^  digital  signal 
processing  technologies,  it  is  possible  to  develop  digital  microwave  receivers  which  can 
meet  this  requirement.  A  simple  in— phase/quadrature  (I/Q)  demodulator  architecture 
which  can  measure  the  amplitude,  phase  and  instantaneous  frequency  of  radar  signals  is 
examined  in  this  report.  The  attractive  features  of  this  architecture  for  radar  ESM 
applications  are:  (i)  wide  instantaneous  bandwidth  because  the  negative  and  positive 
frequencies  with  respect  to  the  local  oscillator  frequency  can  be  distinguished,  and  (ii) 
simple  algorithms  for  the  extraction  of  modulation  characteristics  can  be  used  so  that 
nearly  real-time  results  can  be  obtained.  However,  because  the  splitting  of  the  signal  into 
its  in-phase  and  quadrature  components  is  implemented  using  analog  components,  there 
will  be  amplitude  and  phase  imbalances  between  the  two  channels  and  DC  offsets,  which 
will  introduce  systematic  errors  to  the  measurement.  Two  techniques,  compensation  and 
low— pass  filtering  using  a  moving  average,  are  proposed  in  this  report  for  improving  this 
simple  architecture. 

The  effect  of  amplitude  and  phase  imbalances  and  DC  offsets  is  found  to  produce 
ripples  on  the  measurement  and  the  frequency  of  the  ripples  is  harmonically  related  to  the 
baseband  frequency  of  the  signal.  For  a  given  set  of  imbalance  errors  and  DC  offsets,  the 
error  introduced  by  the  ripples  on  the  instantaneous  frequency  is  more  severe  and  is 
directly  proportional  to  the  baseband  frequency  while  the  error  is  constant  for  both  the 
amplitude  and  phase. 

Some  I/Q  networks  have  been  evaluated  and  the  mean  imbalance  errors  and  DC 
offsets  over  a  large  frequency  range  can  be  large.  In  order  to  keep  the  systematic  errors 
small,  some  form  of  compensation  is  needed  to  remove  the  mean  errors  and  offsets.  Once 
the  mean  errors  and  offsets  are  eliminated,  the  residual  RMS  errors  as  a  function  of 
frequency  and  input  power  level  are  quite  small.  If  further  improvement  is  required, 
calibration  dependent  on  frequency  and  input  power  levd  may  be  needed. 

The  measurement  accuracy  on  the  amplitude,  phase  and  instantaneous  frequency 
has  also  been  analyzed  and  given  in  terms  of  input  SNR.  For  large  input  SNR,  the  standard 
deviation  of  the  phase  error  decreases  inversely  proportional  to  the  square  root  of  the  input 
SNR  while  the  standard  deviation  of  the  envelope  remains  constant  as  the  input  signal 
level  is  varied. 

There  are  two  basic  functions  performed  by  using  a  moving  average.  The  effective 
noise  bandwidth  of  the  I/Q  network  can  be  reduced  with  no  appreciable  degradation  of  the 
signal  of  interest,  which  will  improve  the  output  SNR.  The  ripples  produced  by  the 
imbalance  errors  and  DC  offsets  are  usually  outside  the  video  bandwidth  of  the  signal  of 


V 


interest  and  thus  can  be  reduced  by  low— pass  filtering.  In  addition,  a  moving  average  is 
easy  to  implement  digitally  and  does  not  introduce  distortions  such  as  ringing  and 
overshoot. 

Both  the  techniques  of  compensation  and  low— pass  filtering  using  a  moving 
average  have  been  successfully  demonstrated  on  the  demodulation  of  pulsed  and  linear  FM 
signals.  They  have  been  shown  to  be  effective  for  improving  the  output  signal— to-noise 
ratio  and  reducing  systematic  errors. 
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1.0  INTRODUCTION 


Due  to  the  increasing  density  and  complexity  of  radar  signal  waveforms,  it  is  liO 
longer  easily  possible  to  sort  out  uniquely  and  identify  each  radar  emitter  using 
conventional  signal  parameters  such  as  pulse  width,  radio  frequency  (RF),  amplitude  and 
pulse  repetition  frequency.  As  a  result  there  is  a  requirement  for  a  radar  Electronic  Support 
Measures  (ESM)  receiver  to  measure  the  modulation  characteristics  of  radar  signals  and  to 
provide  additional  parameters  which  can  be  used  to  identify  unambiguously  each  type  of 
radar  emitter.  With  the  advent  of  fast  A/D  converters  and  high— sp^  digital  signal 
processing  technologies,  it  is  possible  to  develop  digital  microwave  receivers  which  can 
meet  this  requirement  [1].  A  simple  in— phase/quadiature  (I/Q)  demodulator  architecture 
which  can  measure  the  amplitude,  phase  and  instantaneous  frequency  of  radar  signals  is 
shown  in  Fig.  1.  The  attractive  features  of  this  architecture  for  radar  ESM  applications  are. 
(i)  wide  instantaneous  bandwidth  because  the  negative  and  positive  frequencies  with 
respect  to  the  local  oscillator  frequency  can  be  distinguished,  and  (ii)  simple  algorithms  for 
the  extraction  of  modulation  characteristics  can  be  used  so  that  nearly  real-time  results 
can  be  obtained.  However,  because  the  splitting  of  the  signal  into  its  in— phase  and 
quadrature  components  is  implemented  using  an^dog  components,  there  will  be  amplitude 
and  phase  imbalances  between  the  two  channds  and  DC  offsets,  which  will  introduce 
systematic  errors  to  the  measurement [2—3]. 

This  report  addresses  this  I/Q  sampling  architecture  for  the  measurement  of 
amplitude,  phase  and  instantaneous  frequency  of  radar  signals.  A  time-domain  analysis,  on 
the  effects  of  phase  and  amplitude  imbalances  between  the  I  and  Q  channels  and  DC 
offsets,  is  given.  A  calibration  technique  which  can  be  used  to  compensate  for  the 
imbalances  and  offsets  is  discussed.  The  effect  of  noise  on  the  measurement  is  analyzed. 
Imbalance  errors  and  DC  offsets  of  I/Q  networks  are  measured  and  analyzed.  Finally,  a 
post— processing  technique  employing  moving  averages,  which  can  be  us^  for  improving 
the  output  signal— to— noise  ratio  and  reducing  systematic  errors,  is  also  presented. 


2.0  I/Q  DEMODULATION 

An  incoming  RF  signal  is  usually  down— converted  to  an  intermediate  signal  (IF) 
before  it  is  applied  to  the  I/Q  network.  The  IF  signal  to  the  I/Q  demodulator,  as  shown  in 
Fig.  1  ,  can  be  expressed  in  the  form  of 

8(t)  =  a(t)cos[4;^t  +  ^(t)]  (1) 

where  a(t)  is  the  amplitude  or  envelope,  is  the  angular  carrier  frequency  and  |^(t)  is  the 

phase  function  of  the  signal.  It  is  first  bandpass  filtered  and  then  equally  power-divided 
into  two  paths.  The  signal  in  each  path  is  then  mixed  down  to  baseband  by  the  use  of  a 
local  osciUator  signal.  The  two  local  oscillator  signals  are  derived  from  the  same  source,  but 
are  90  degrees  out  of  phase.  The  resultant  in-phase  and  quadrature  baseband  signals  after 
the  low— pass  filter  (LPF)  are. 


Si(t)  =  K/2  a(t)cos[  6>io)t  +  ^(t)  ~  7]  =  K/2  a(t)cos[^(t)] 


(2) 


and 

Sq(t)  =  K/2  a(t)bin[^(t)]  (3) 

respectivdy,  where  K  is  the  net  gain  in  each  path  ,  wio  is  the  angular  frequency  of  the  local 
oscillator  with  initial  phase  7,  /?(t)  is  the  phase  function  of  the  baseband  signal  and  a 
constant  delay  introduced  in  each  path  has  been  neglected.  In  this  ideal  case,  the  two 
channels  have  been  assumed  to  be  perfectly  matched  in  amplitude,  90  degrees  out  of  phase 
and  with  no  DC  offsets. 

The  instantaneous  power  of  the  envelope  of  the  input  signal  is  simply  related  to  its 
in— phase  and  quadrature  baseband  components  by 

a»(t)  =  4/K“  (  s?(t)  +  sj(t)l  ,  (4) 

the  signal  phase  function  is  given  by 

^(t)  =  ^(t)  -  (Wg-  4;io)t  +  7 

=  tan*‘  [Sq(t)/Si(t)]  -  (Wg-  wio)t  +  7  (5) 

and  the  instantaneous  angular  frequency  is 

2ri(t)=  ^  ^  [  tan-  lS,(t)/Si(t)l  ]  -  («.-  «.).  (6) 

The  in-phase  and  quadrature  baseband  signals  are  usually  sampled  and  digitized  at  tn  = 
nTs  +  Tg,  where  Ts  is  the  sampling  interval,  T^  is  the  initial  time  and  n  =  0,1,2,...  .  In 

this  case,  the  sampled  instantaneous  angular  frequency  is  then  approximately  given  by 

2rf(t„)»  [^(l„)-^(t„.,)]/T,  (7) 

3.0  AMPLITUDE  AND  PHASE  IMBALANCES  AND  DC  OFFSETS 

In  the  implementation  of  the  I/Q  network,  there  will  be  differential  gain,  DC  offsets 
and  phase  deviation  from  the  ideal  90  degrees  between  the  two  channels  [3].  Equations  (2) 
and  (3)  can  then  be  rewritten  in  a  more  general  form  as 

Si(t)  =  Ki/2  a(t)cos[/J(t)  +  ]  +  aio  (8) 

and 

Sq(t)  =  Kq/2  a(t)sin[^(t)  +  ^q  ]  +  aqo  (9) 

respectively,  where  aio  and  aqo  are  the  amplitude  DC  offsets,  Ki  and  Kq  are  the  gains,  and 
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and  4q  &re  the  phases  of  the  in— phase  and  quadrature  channels  respectively,  aio,  aqo,  Ki, 
Kq  and  are  in  general  a  function  of  frequency.  For  narrow— band  signals,  they  can  be 
assumed  to  be  approximately  constant. 

3.1  Envelope  Measurement 

Substituting  Eqs.(8)  ai  .d  (9)  into  Eq.(4),  the  instantaneous  power  of  the  envelope  of 
the  baseband  sign^  can  be  shown  to  be  [4] 

Sh)  +  Sj(t)  =  [a(t)Ki]V4  {  1/2  +  {aio/[a(t)Ki/2]}2+  RV2  +  {aqo/[a(t)Ki/2]}2 
+  ^  cos[^(t)  +  ]{2  aio/[a(t)Ki/2]  +  2  R  aqo/[a(t)Ki/2]  sin[(^q  -  )]} 

+  2  R  aqo/[a(t)Ki/2]  sin[/?(t)  +  ^i]  cos[(^q  -  )]  j 

+  1/2  cos  {  2  [^(t)  +  ]}  -RV2  [  cos  {  2  +  ^i]  +  [  2  (  ^q-  )]  }  ]| 

(10) 


The  envelope  has  been  expressed  in  terms  of  normalized  parameters,i.e., 
fi)  amplitude  imbalance  ratio  fR  =  Kq/Ki),  (ii)  phase  imbalance  (^jj  —  ),  and 

(iii)  normalized  DC  offsets  aio/[a(t)Ki/2]  and  aqo/[a(t)Ki/2]  of  the  in-phase  and 
quadrature  channels  respectively.  If  the  envelope  oi  the  input  signal  [a(t)]  is  constant, 
then  there  are  basically  three  groups  of  terms  inside  the  braces;  DC  terms,  fundamental 
baseband  frequency  terms  and  second  harmonic  baseband  frequency  terms.  The  DC  terms 
is  a  function  of  the  amplitude  imbalance  and  DC  offsets,  independent  of  phase  imbalance. 
A  DC  offset  wiU  introduce  a  fundamental  harmonic  frequency  which  has  the  same 
frequency  as  the  baseband  signal.  Other  errors  wiU  only  affect  the  magnitude  of  these 
ripples.  Without  DC  offrets,  ripples  of  the  fundamental  frequency  will  disappear 
completely.  Ripples  with  a  second  harmonic  frequency  will  appear  only  when  there  is  an 
amplitude  or  phase  imbalance.  When  the  two  channels  are  matched,  all  the  ripples  will 
disappear  and  the  terms  inside  the  bracdets  will  be  unity.  The  effects  of  the  three  different 
sources  of  errors  on  the  measurement  have  been  analyzed  in  detail  [4].  The  ratios  of  the 
peak— to— peak  ripple  to  the  mean  are  plotted  in  Figs.2  and  3,  for  tne  cases  when  there  are 
only  amplitude  and  phase  imbalances,  and  when  there  are  only  DC  offsets  respectively. 
When  there  are  only  amplitude  and  phase  imbalances  in  the  network,  the  ratio  of 
peak— to-peak  ripple  to  mean  is  a  function  of  the  absolute  phase  imbalance  and  the 
amplitude  imbalance  ratio  expressed  in  dB.  In  the  case,  when  there  are  only  DC  offset 
errors,  the  ratio  is  just  a  function  of  the  absolute  value  of  the  normalized  DC  offsets. 


3.2  Phase  and  Frequency  Measurement 

Substituting  Eqs.(8)  and  (9)  into  Eq.(5),  results  in  the  phase  of  the  baseband  signal 
given  by 

o(t)  =  tan'*  I  R  1^  sin[^(t)  +  ^i]  cos((^q  -  )]  +  cos[^(t)  +  ^i] 

•  8in[(^q  -  )]  +  aqo/[a(t)Kq/2lj/  ^  cos[;3(t)  +  ]  +  aio/[a(t)Ki  /2]j  |  (11) 
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Fig.  3  Normalized  Peak-to— Peak  Ripple  Versus  DC  Offsets 
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where  the  argument  has  also  been  expressed  in  terms  of  normalized  parameters. 

The  deviation  in  measurement  from  the  case  of  an  ideal  sampling  network  is 
emphasized  in  this  analysis.  The  phase  error  is  defined  bv  the  difference  between  the 
measured  phase  function  and  the  input  phase  function  or  the  baseband  signal  as 

Aa(t)=  o(t)-ajt)  (12) 

where 

flo(‘)  =  ^(0  +  (13) 

is  the  input  phase  function  of  the  baseband  signal  with  a  constant  phase  ofiset  i^i.  For  a 
specific  set  of  imbalance  errors  and  DC  o&ets,  the  phase  error  can  be  fully  characterized 
by  plotting  Aaft)  over  an  input  phase  change  of  2t  radians  [4].  The  measured  phase  value 
[a(tj]  is  then  ootained  from  Eq.(12)  once  the  phase  function  of  the  input  signal  is  known.  It 
IS  to  be  noted  that  the  imbalance  errors  and  DC  offsets  are  usually  a  function  of  frequency. 
The  measured  frequency  deviation  from  the  input  baseband  signal  is  then  obtained  by 
differentiating  Eq.(12)  with  respect  to  time 


5[Aa(t)]  5[Aa(t)]  5[a^(t)] 

- ""  51a";(t)]  ”  di  ~  (1^) 

where  5[Aa(t)]/d[ag(t)]  is  the  partial  derivative  of  the  measured  phase  change  with  respect 

to  the  input  phase  and  is  dimensionless.  It  is  also  a  very  useful  parameter  in  characterizing 
the  normalized  frequency  deviation  due  to  the  imbalance  errors  and  DC  offsets.  It  is  found 
that  the  imbalance  errors  and  DC  offsets  have  a  similar  effect  on  both  the  measured  phase 
and  instantaneous  frequency  as  in  the  case  of  the  measured  envelope  [4].  The  major 
difference  is  that  the  ripples  generated  are  more  complex  in  shape.  Figure  4  shows  the 
normalized  maximum  positive  and  maximum  negative  frequency  deviations  as  a  function  of 
the  amplitude  and  phase  imbalances  while  Fig.  5  shows  the  frequency  deviation  due  to  DC 
offsets  only.  In  Fig.  4,  the  maximum  normalized  frequency  deviation  is  only  a  function  of 
the  absolute  phase  imbalance  and  similarly  in  Fig.  5,  it  is  only  a  function  of  the  absolute 
value  of  the  DC  offset. 

For  a  given  set  of  imbalance  errors  and  DC  offsets,  the  normalized  peak  frequency 
deviations  are  completely  determined  by  a  plot  of  d[Aa(t)]/d[a^(t)]  over  an  input  phase 

change  of  2r  radians.  The  measured  maximum  frequency  deviations  from  the  input  signal  is 
determined  by  multiplying  the  normalized  peak  frequency  deviations  by  the  instantaneous 
frequency  of  the  input  signal  at  the  baseband  frequency  as  given  in  Eq.(14).  If  the 
imbalance  errors  and  DC  offsets  are  about  the  same  as  a  function  of  IF  frequency,  the 
measured  frequency  error  due  to  the  ripples  can  be  minimized  by  choosing  an  IF  signal 
which  gives  the  lowest  baseband  frequency. 
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Fig.  5  Maximum  Instantaneous  Frequency  Deviation  Versus  DC  Offsets 


9 


4.0  COMPENSATION  THROUGH  CALIBRATION 

Since  most  of  the  imbalance  errors  and  DC  offsets  are  stationary  and  are  only  a 
function  of  frequency,  some  compensation  can  be  carried  out  to  eliminate  them.  When  a 
calibration  signal  of  the  form  mven  by  Eq.(l)  is  applied  to  the  I/Q  network,  the  resultant 
in— phase  and  quadrature  baseband  signals  are  given  by  Eqs.  (8)  and  (9)  respectively.  When 
the  sine  function  of  the  baseband  signal  in  the  quadrature  channel  is  expressed  in  terms  of 
a  cosine  function  and  the  phase  imbalance,  Eq.(9)  can  be  rewritten  as 

Sq(0  “  ^q/2  a(t)sin[^(t)  +  ^q]  +  aqo 

=  a(t)  Kq/2  ^cos  [/?(t)  +  -  r/4]  cos[(^q  -  ^i)] 

+  cos  [/3(t)  +  ^i]  8in[  (^q  -  ^i)]j  +  aqo  (15) 

Taking  the  Fourier  transform  of  the  above  expression  over  a  finite  duration  T,  we 

have 

.?[Sq(t)]  =  Kq/2  cos[(^q-  ^i)]  exp(-jT/4)  Jf[^a(t)cos  [^(t)  +  ^i]j 
4-  Kq/2  sin[(^)q  -  ^i)]  J>^[^a(t)cos  [^(t)  +  ^i]j  +  .5^[aqo] 

=  -  jKq/2  exp[j(^q  -  ^i)]  3^^  a(t)cos  [/3(t)  +  ^i]j  +  ^[aqo] 


(16) 


T/2 

where  ^[h(t)]  =  J"  h(t)  exp(— j2Tft)  dt  (17) 

-T/2 

is  the  finite  Fourier  transform.  For  a  given  baseband  signal,  the  time  duration  T  must  be 
chosen  so  that  there  will  be  negligible  overlapping  between  the  frequency  spectrum  of  the 
signal  and  the  spectrum  of  the  truncated  DC  o&et. 

The  Fourier  transform  of  the  in— phase  component  is 

.5f[Si(t)]  =  Ki/2  a(t)  cos[^(t)  +  ]]  +  ^[aio]  (18) 

The  DC  frequency  component  due  to  the  DC  offset  can  be  separated  and  removed  if 
there  is  negligible  overlapping  between  the  frequency  components  generated  by  the 
truncated  DC  offset  and  those  due  to  the  signal.  The  DC  offsets  of  the  in-phase  and 
quadrature  channels  are  simply  related  to  its  finite  Fourier  transform  by 

aio  =  ^  <?’[aio]  at  f  =  o|  /T  (19) 
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and 


V  =  [  <^[^0]  at  f  =  oj  /T  (20) 

respectively.  Taking  the  ratio  of  the  Fourier  transform  of  the  simal  in  the  quadrature 
channel  to  the  Fourier  transform  of  the  signal  in  the  in— phase  channel,  we  have 

^[Sq(t)]/ .5^[Si(t)]  =  -  j  exp[j(^q  -  ^i)]  Kq/Ki 


=  R  exp  ^j((^q  -  ^i)  -  t/4]  j  (21) 

Once  the  ratio  of  «?’[Sq(t)J/ J^[Si(t)]  is  computed,  the  relative  amplitude  imbalance  ratio  R 
and  the  phase  imbalance  (^q  —  as  a  function  of  frequency  can  then  be  determined.  Using 
the  above  procedure  and  knowing  the  input  IF  frequency,  all  the  imbalance  errors  and  DC 
offsets  of  the  network  can  be  measured. 

In  practice,  the  baseband  signals  are  sampled  and  digitized.  As  a  result,  a  discrete 
Fourier  transform  (DFT)  has  to  be  used  and  t  is  replaced  by  tn-  Moreover,  a  CW 
calibration  signal  is  used  for  two  main  reasons:  (i)  precise  CW  sources  are  readily  available 
and  (ii)  the  number  of  data  samples  generated  at  a  particular  frequency  of  interest  can  be 
arbitrarily  large.  A  large  number  of  data  samples  for  the  DFT  can  improve  the  output 
signal— to— noise  ratio  and  thus  the  ultimate  accuracy  of  the  measurement.  Since  both  the 
input  IF  frequency  and  the  resultant  baseband  frequency  are  known,  the  number  of  sample 
points  can  also  be  chosen  to  satisfy  the  criterion  of  having  a  truncation  interval  for  the 
DFT  equal  to  an  integer  multiple  cd  the  period  of  the  baseband  signal.  Once  this  criterion 
is  satisfied,  the  leakage  problem  due  to  time  domain  tnmcation  associated  with  DFT  can 
be  eliminated  [5].  This  is  due  to  the  fact  that  when  the  truncation  interval  is  equal  to  an 
integer  multiple  of  the  period  of  the  baseband  signal,  the  CW  signal  will  locate  at  one  of 
the  nulls  of  the  Fourier  transformed  sine  function.  As  a  result,  no  distortions  axe 
introduced  to  either  the  CW  signal  or  the  DC  offset.  The  DC  offset  is  determined  by 
dividing  the  DC  component  of  the  DFT  by  the  total  number  of  sample  points  used. 

Once  the  imbalance  errors  and  DC  ofisets  as  a  function  of  frequency  are  known  the 
errors  can  be  eliminated  by  a  calibration  technique.  In  the  simplest  case  where  the 
imbalance  errors  and  DC  offsets  of  the  I/Q  network  over  the  signal  bandwidth  are 
approximately  constant,  a  single  set  of  (Vibration  values  can  be  applied  to  every  data 
point  once  the  approximate  frequency  of  the  signal  is  estimated. 

Recall  that  the  objective  is  to  obtain  a^(t)  in  Eq.(13).  For  this  we  need  expressions 

for  a(t)8in  [/J(t)  +  and  a(t)  cos[/J(t)  +  ^i].  Using  both  Eqs.(8l  and  (9),  a(t)sin  [^(t)  +  ^i] 
can  be  expressed  in  terms  of  the  measured  baseband  signals,  ana  the  imoalance  errors  and 
DC  offsets  as 

a(t)sin  [/?(t)  +  ^i]  =  ^[Sq(t)  -  aqoJ/(  Kq/2  )  -  [S.(t)  -  aio]/(Ki/2  ) 

•  8in[  (^q  -  ^i)lj  /cos[  (^q  -  ^i)]  (22) 
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Similarly,  from  Eq.(18),  we  also  have 

a(t)  cos[/9(t)  +  ^i]  =  [S.(t)  -  aio]/(Ki/2)  (23) 

Dividing  Eq.(22l  by  Eq.(23)  and  taking  the  arctangent  of  the  ratio,  the  phase  function  of 
the  compensated  baseband  signal  is  given  by 

[/?(*)  +  ^i]  =  tto(t) 

■  ^ 

=  tan-‘  [^[Sq(t)  -  aqo]/R  -  [S.(t)  -  aio]sin{(^q  -  ^i)]]  /  |^cos[(^q  -  ^i)]  [S.(t)  -  aio]] 

■  ^ 

(24) 

Since  the  compensated  phase  function  of  the  baseband  signal  is  known,  the 
compensated  amplitude  of  the  envelope  is  obtained  by  using  Eq.(23).  Hence 

a(t)  =  2/Ki  [S.(t)  -  aio]/cos[^(t)  +  ^i]  (25) 

A  functional  block  diagram  for  the  above  compensation  and  computation  process  is 
shown  in  Fig.  6. 

The  next  section  considers  the  effect  of  input  noise  on  the  measurements. 

5.0  SIGNAL  PLUS  NOISE 

In  radar  ESM  applications,  the  intercepted  radar  signal  is  usually  amplified  by  a 
pre-amplifier  at  the  RF  and  then  down-converted  to  an  IF  signal.  The  IF  signal  is  then 
amplified  again  before  it  is  applied  to  the  input  of  the  I/Q  network.  The  input  noise 
contribution  is  usually  dominated  by  the  RF  front-end  pre— amplifier  noise  since 
subsequent  stages  in  the  receiver  are  preceded  by  sufficient  gain  in  order  to  overcome 
additional  noise.  The  active  components  in  the  I/Q  network  will  also  introduce  noise.  In 
addition,  quantization  noise  will  also  be  introduced  when  the  baseband  signals  are 
digitized.  It  is  assumed  that  the  RF  noise  is  dominant  and  is  white  Gaussian  with  zero 
mean  and  spectral  density  if (2  at  the  input  of  the  I/Q  network.  After  mixing  down  to 
baseband  and  with  low  pass  filtering,  the  noise  in  tne  two  channels  are  uncorrelated  Avith 

each  having  a  spectral  density  K’>//4  .  For  Bjj,  =  2  Bj^p,  where  Bj^p  is  the  LPF 

bandwidth,  and  if  the  signal  is  approximately  located  at  the  center  of  the  band  pass  filter 
(BPF),  the  mean  si^al— to-noise  ratio  (SNR)  at  the  output  of  each  channel  is  about  the 
same  as  that  at  the  input.  The  signals  in  both  the  in— phase  and  quadrature  channels  are  of 

equal  peak  amplitude  {Ka(t)/2].  The  mean  SNR  [a^(t)/(2(r2)]  channel  is  equal  to 

5.1  Phase  Error  As  a  Function  of  Signal— to— Noise  Ratio 

At  time  t  =  t^,  when  the  phase  of  the  input  signal  0^  =  /?(t^),  the  probability 


density  function  as  a  function  of  SNR  can  be  shown  to  be  [4,6] 

U  (^)=  l/(2T)exp[-SNR] 

+  (SNR)*/’  cos(^  -  ^^)/[2(T)‘/’].exp  [-  SNR  [sin{9  -  Jj’j 
.  [l  +  erf[(SNR) */’  cos( -  ^^)]  j  ,  for  ]  ]  <  x  (26) 

where 

X 

erf(x)  =  2/Vt  J*  exp(—  t’)  dt  (27) 

0 

For  small  SNR  ratios,  the  second  term  is  negligible  and  the  probability  function  is 
approximately  given  by 

(^)  ®  1/(2t)  I  for  l^l<  (28) 

which  is  uniform  in  distribution.  For  large  SNR  ratios,  the  first  term  is  negligible  and 
erf[(SNR)‘/’cos(tf-^^)]a  1 

co8{0  -  tf^)  a  1  and 

sin{9-9J»  (29) 


Substituting  the  above  approximations  into  Eq.(26),  we  have 


fg  (0)  a  (SNR)*/ V(x)*/’  exp  [-  SNR  (0  -  0J^  ,  for  ]  ^  1  <  x 
0 

or  »  1  /((2r)'/ Va(‘o)I  «p  [-  («  -  «,)V{2('’/»(tJl'}]  ,  for  I « I  <  r 


(30) 

which  is  approximately  Gaussian. 

Once  f^  (0)  is  determined,  the  mean  and  variance  of  the  measured  phase  angle  can 

O 

easily  be  determined  by  using  the  following  equations: 


E  (»,)  =  Vf 

o 


QD 


0  f^  (0)  d0,  where  |  ^1<  x 

O 


(31) 
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and 


Variance  = 


0 


where  \  9\<t 


(32) 


Substituting  Eq.(26)  into  £qs.(31)  and  (32),  the  mean  and  variance  of  the  phase  can 
be  estimated.  For  small  signal-to-noise  ratios,  the  mean  {ri^  )  is  0  and  the  standard 

O 

deviation  {a ^  )  is  radian. 


For  large  signal— to-noise  ratios,  the  mean  phase  value  is  simply  Oq  and  the 
standard  deviation  is  tf’/a(t^)  or  1/(2  SNR)*/*  radian. 


5.2  Amplitude  Error  As  a  Function  of  Signal— to— Noise  Ratio 

At  time  t  =  t^,  when  the  amplitude  of  the  input  signal  a(t)  =  a(t^),  the  probability 
density  function  of  the  envelope  is  given  by[4,6] 

f.(z')  =  *7<'«p[-(2'’  +  2SNR)/2j  y  (2  SNR)*/’ 2']  ,  for  2'  >0 

(33) 

where  <r  is  the  standard  deviation  of  the  noise  in  either  the  I  or  Q  channel  and  z'  =  z/tr  is 
the  normalized  amplitude  of  the  envelope  and  is  the  modified  Bessel  function  of  order 

zero  defined  by 

2t 

IoW  =  ^/(2^)  /  exp[xcos(^)]  d^  (34) 

0 

The  mean  of  the  envelope  is  given  by 


E  (z')  =  j/z' 


where  z'  >  0 


(35) 


and  the  variance  is 
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(36) 


=  j  {z'  —  itz'y  dz\  where  z'  >  0 

—  OD 

For  small  signal— to-noise  ratios, 

{z(z')  a  exp(—z'^/2)  ,  for  z' >  0  (37) 

which  is  of  Rayleigh  distribution  and  the  mean  (tr  i/t')  is  and  the  standard 

deviation  (itz)  is  (2  —  rf  2y^^  a. 

For  large  signal— to-noise  ratios, 

I^[(2  SNR)'/’  z']  s  exp  [  (2  SNR)'/’  z']/[2t  (2  SNR)'/’  z')'/’ 


(38) 


Substituting  Eq.(38)  into  Eq.(33),  the  probability  density  function  is  approximately 
given  by 


fa(z')  a  z'/{(r  [25r  (2  SNR)‘/^  s^pj^—  [z'  —  (2  SNR)*^^  /2j  ,  for  z'  >  0 


or 

a  z7[(r[2ra(t,)/irz']‘/’|  expj^- [z' - a(t^)/x  ]*/2j  , 


(39) 


for  z'  >  0 

Equation  (39)  is  approximately  a  Gaussian  distribution  with  mean  equal  to  a.{t^)/<r. 
At  z'  a  z,{t^)l(r,  the  coefficient  of  the  above  expression  is  approximately  equal  to 

l/[(2T)*/’fl-]  and  thus  the  standard  deviation  of  the  envelope  is  <r. 

6.0  MOVING  AVERAGE 

A  moving  average  is  performed  on  the  processed  data  of  the  envelope,  phase  and 
instantaneous  frequency  of  the  signal  in  the  time  domain.  This  is  carried  out  by  adding  N 
contiguous  samples  and  an  average  is  taken  on  the  accumulated  sum.  The  same  operation 
is  repeated  on  the  next  N  samples.  The  transfer  function  is  given  by[7] 


H(f)/N  =  sin(Nrf  Tg)/  [N  sin(xf  T^)]  exp  [  -j(N  -  1)  xf  T  J  (40) 
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The  absolute  value  of  the  transfer  function  is  plotted  in  Fig.  7  for  various  values  of 
N.  As  can  be  seen  from  the  plot,  the  operation  of  moving  average  on  the  sampled, 
processed  data  has  the  same  effect  as  low— pass  filtering  in  the  fiequency  domain.  There  are 
two  basic  functions  performed  by  using  the  moving  average.  The  first  function  is  to  reduce 
the  video  bandwidth  and  thus  the  effective  noise  bandwidth  [8]  of  the  I/Q  network  is 
reduced  with  a  corresponding  improvement  on  the  output  SNR.  The  second  function  is 
when  there  are  amplitude  and  phase  imbalances  and  DC  offsets  in  the  I/Q  network  , there 
will  be  ripples  introduced  on  the  amplitude,  phase  and  instantaneous  frequency  as  pointed 
out  in  Section  3.0.  Depending  on  the  frequency  of  the  baseband  signal,  these  ripples  can  be 
outside  the  video  bandwidth  of  the  signal  of  interest  and  thus  reduced  by  low— pass 
filtering.  The  moving  average  is  easy  to  implement  digitally  and  does  not  introduce 
distortions  such  as  ringing  and  overshoot  as  other  types  of  filters  [9].  As  given  in  Fig.  7, 
the  sidelobe  levels  are  determined  by  Eq.(40)  as  a  function  of  N.  For  large  N,  the  first  peak 
sidelobe  level  is  13.46  dB  down  and  the  peak  sidelobe  level  decreases  directly  proportional 
to  1/f.  If  it  is  desirable  to  have  lower  sidelobe  levels,  a  multi-layer  moving  average  can  be 
performed.  Fig.8  shows  the  resultant  effect  when  two  layers  are  used.  In  tUs  case  the  data 
is  filtered  by  a  moving  average  of  N  =  7  and  followed  by  N  =  5. 


7.0  EXPERIMENTAL  MEASUREMENTS 

Implementation  of  the  I/Q  network  is  carried  out  by  the  use  of  a  quadrature  mix'^, 
a  local  oscillator  and  a  dual— channel  digital  scope.  The  quadrature  mixers  used  are  from 
Anaren  Microwave  Inc.  and  KDI/Triangle  Electronics  Inc.  The  digital  scope  is  either  a 
10— bit  or  8— bit  dual— chaimel  scope  from  LeCroy  Corporation.  On  any  given  set  of 
measurements,  the  set-up  is  usua^y  ”  warmed  up  ”  for  at  least  20  minutes  before  the 
testing  begins.  Depending  on  the  type  of  testing,  the  actual  time  taken  to  complete  the  test 
can  last  over  a  time  period  of  tens  of  minutes  and  sometimes  even  hours.  As  a  result,  if  the 
characteristics  of  the  I/Q  network  varies  as  a  function  of  time,  they  are  alsr  reflected  in  the 
experimental  data. 

7.1  Imbalance  Enors  and  DC  Offeets  of  I/Q  Network 

The  imbalance  errors  and  DC  offsets  of  I/Q  network  are  measured  by  using  the 
technique  described  in  Section  4.  The  I/Q  network  includes  a  quadrature  mixer,  a  3-dB 
pad  ,  an  intercoimecting  cable  and  the  10— bit  LeCroy  9430  duzd-channel  digital  scope.  The 
3— dB  pad  is  used  to  improve  the  impedance  matching  between  the  quadrature  mixer  and 
the  scope.  The  in— phase  and  quadrature  signals  are  sampled  and  di^tized  by  the  scope  at  a 
100— MHz  sampling  rate  with  a  low— pass  filter  bandwidth  of  approximately  150  MHz. 
DiHerent  quadrature  mixers  are  tested  as  a  function  of  the  LO  Requency  ,  LO  power  level 
and  input  signal  power  level.  A  CW  signal  from  a  frequency  synthesizer  is  applied  to  the 
input  of  the  I/Q  network  and  is  then  adjusted  to  produce  a  baseband  signal  from  —  50  MHz 
to  50  MHz  at  1— MHz  steps.  A  4000-point  DFT  is  performed  on  the  digitized  baseband 
signals  to  extract  the  relative  amplitude  and  phase  diHerence  between  the  two  channels  and 
the  DC  offset  in  each  channel.  A  typical  signal— to— noise  ratio  of  approximately  48  dB  is 
achieved  on  each  sample.  When  a  DFT  is  performed  on  4000  sample  points,  the 
signal— to— noise  ratio  is  approximately  improved  by  an  additional  factor  of  36  dB.  As  a 
result,  the  combined  signal— to— noise  ratio  is  approximately  84  dB.  The  noise  distribution 
is  approximately  Gaussian  and  the  errors  introduced  by  the  noise  as  discussed  in 
Section  5.0  are  negligibly  small  and  its  effect  on  the  measurements  can  be  neglected. 


The  amplitude  ratio  R  expressed  by  20*logio(Kq/Ki)  in  dB  versus  baseband 
frequency  is  plotted  in  Fig.9  (a)  for  a  typical  I/Q  network.  The  quadrature  mixer  used  is 
from  Anaren  Microwave  Inc.  The  LO  frequency  is  set  at  1.5  GHz  with  a  power  output  of  10 
dBm.  The  input  signal  power  level  at  1.5  GHz  is  adjusted  so  that  approximately  80%  of 
full  scale  deflection  is  achieved  on  the  digital  scope,  which  is  equivalent  to  ±  32  mV  in  a 
50-fl  termination.  The  same  input  power  level  is  maintained  as  the  signal  frequency  is 
varied.  The  mean  amplitude  ratio  over  the  full  baseband  frequency  range  of  ^  50  MHz  is 
computed  to  be  1.034  (0.291  dB).  The  RMS  value  is  0.011  and  when  expressed  in  terms  of 
20*logio  [(  mean  +  RMS)/(mean  —  RMS)],  it  is  0.186  dB. 

The  absolute  phase  difrerence  between  the  quadrature  channel  and  the  in— phase 

channel  is  plotted  in  Fig.9  (b).  The  mean  value  is  90.46°  and  the  RMS  value  is  1.81°.  As 
can  be  seen  from  the  plot,  the  phase  difference  varies  approximately  linearly  as  a  function 
of  baseband  frequency.  As  a  result,  the  linear  component  can  be  removed  by  adjusting  the 
cable  lengths  in  both  the  in— phase  and  quadrature  channels.  By  applying  a  least  square  fit 
to  the  data  in  Fig.  9(b),  the  absolute  phase  difference  at  zero  baseband  frequency  is  found 

to  be  90.465°,  the  slope  is  —  0.0606°/MHz  and  the  resultant  RMS  value  is  0.506°  over  the 
baseband  frequency  range  of  ±  50  MHz.  If  the  frequency  range  is  reduced  to  ±  30  MHz,  the 

absolute  phase  difference  at  zero  Hz  is  90.365°,  slope  is  —  0.08075°/MHz  and  the  RMS 

value  is  reduced  to  0.262°. 

The  DC  offsets  in  mV  of  both  the  in-phase  and  quadrature  channels  are  plotted  in 
Fig.9  (c).  The  mean  DC  offset  values  of  the  in-phase  and  quadrature  channels  are  1.78  mV 
and  —  2.29  mV  respectively.  The  RMS  values  for  both  the  in— phase  and  quadrature 
channels  are  0.064  mV  and  0.073  mV  respectivdy  over  the  baseband  frequency  range  ■'f 
±  50  MHz. 

The  measured  power  levels  in  both  the  in-phase  and  quadrature  channels  are 
plotted  in  Fig.9(d)  versus  baseband  frequency. 

Using  the  same  set-up  but  with  the  input  signal  power  reduced  by  approximately 
4  dB  as  shown  in  Fig.lO  (dl,  the  above  set  of  measurements  is  repeated.  The  purpose  of 
this  measurement  is  to  evaluate  changes  in  the  imbalance  errors  and  DC  offsets  of  the  I/Q 
network  as  the  input  signal  power  level  is  reduced.  The  imbalance  errors  and  DC  offsets 
are  plotted  in  Fig.  10(a)  to  Fig.  10(c).  By  comparing  the  two  sets  of  measurements  over  the 
full  baseband  frequency  range  of  ^  50  MHz,  a  RMS  value  on  the  difference  in  amplitude 

ratio  is  found  to  be  0.00318,  a  phase  difference  of  0.126°  and  DC  o&ets  of  0.137  mV  and 
0.223  mV  in  the  in— phase  and  quadrature  channels  respectively.  As  a  result,  there  are 
ne^gible  differences  in  the  imbalance  errors  and  DC  offsets  when  the  input  power  level  is 
reduced  by  4  dB. 

Another  set  of  measurements  is  taken  to  evaluate  the  imbalance  errors  and  DC 
offsets  over  a  much  wider  input  signal  power  range  for  one  specific  signal  frequency.  The 
signal  frequency  is  set  at  1040  MHz  while  the  LO  frequency  is  at  1  GHz  with  a  power 
output  of  approximately  13  dBm.  The  input  signal  power  level  is  varied  in  3— dB  steps  and 
at  each  step,  the  imbalance  errors  and  DC  ofcets  of  the  I/Q  network  are  computed.  They 
are  plotted  in  Fig.  11(a)  to  Fig.  11(c).  As  the  input  power  is  reduced,  the  output 
signal— to-noise  ratio  is  also  reduced.  In  order  to  compensate  for  the  reduced  output 
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Fig.9  Imbalance  Errors  and  DC  Offsets  of  I/Q  Network 
( Input  Signal  Power  Level  »  —  19  dBm) 
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(b)  Absolute  Phase  Difference  Versus  Baseband  Frequency 


Fig.9  Imbalance  Errors  and  DC  Offsets  of  I/Q  Network 
( Input  Signal  Power  Level  »  —  19  dBm) 
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Fig.  9  Imbalance  Errors  and  DC  Offsets  of  I/Q  Network 
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Fig.  10  Imbalance  Errors  and  DC  Offsets  of  I/Q  Network 
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Fig.  10  Imbalance  Errors  and  DC  Offsets  of  I/Q  Network 
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(a)  Amplitude  Ratio  Versus  Input  Signal  Level 

Fig.  11  Imbalance  Errors  and  DC  Offsets  of  I/Q  Network 
Versus  Input  Signal  Level 
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Fig.  11  Imbalance  Errors  and  DC  Offsets  of  I/Q  Network 
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Fig.  11  Imbalance  Errors  and  DC  Offsets  of  I/Q  Network 
Versus  Input  Signal  Level 


signal— to— noise  ratio,  more  sample  points  are  needed  for  the  DFT.  For  this  set  of 
measurements,  the  number  of  sample  points  has  been  increased  from  4,000  to  40,000.  As 
can  be  seen  from  the  graphs,  the  imbalance  errors  and  DC  offsets  do  vary  slightly  over  the 
24— dB  testing  range. 

Under  different  operating  conditions,  the  imbalance  errors  and  DC  offsets  of  the  I/Q 
network  have  also  been  evaluate.  It  is  found  that  the  imbalance  errors  and  DC  offrets  are 
a  strong  function  of  the  LO  frequency  and  power  level  to  the  quadrature  mixer.  Once  the 
optimum  operating  LO  frequency  and  power  level  are  chosen  and  fixed,  the  imbalance 
errors  and  DC  offsets  are  rdatively  insensitive  to  input  signal  power  level  and  frequency. 

In  summary,  the  mean  imbalance  errors  and  DC  offsets  of  typical  I/Q  networks  over 
a  large  frequency  range  can  be  large.  In  order  to  keep  the  systematic  errors  small,  some 
form  of  calibration  is  needed  to  remove  the  means  of  the  imbalance  errors  and  DC  offsets. 
Once  the  means  are  removed,  the  RMS  values  of  the  imbalance  errors  and  DC  offrets  as  a 
fimction  of  signal  frequency  and  power  is  quite  small.  If  further  reduction  in  the  systematic 
errors  is  required,  a  more  daborate  calibration  is  required  as  a  function  of  the  signal 
frequency  and  may  even  the  signal  power. 

7.2  Amplitude,  Phase  and  Instantaneous  Frequency  Measurement  Versus  Input  SNR 

For  this  set  of  measurements,  a  CW  signal  from  a  frequency  synthesizer  is  amplified 
by  two  high  gain  IF  amplifiers  before  it  is  applied  to  the  I/Q  network.  The  total  noise 
appearing  on  the  scope  is  dc  linated  by  the  noise  from  the  first  amplifier  and  is  of  Gaussian 
distribution  with  zero  mean.  The  signd  and  noise  in  the  in-phase  and  quadrature  channels 
are  sampled  and  digitized  at  a  rate  of  100  MHz.  The  frequency  of  the  CW  signal  from  the 
synthesizer  is  findy  tuned  so  that  the  resultant  baseband  signal  is  almost  completdy 
phase-locked  to  the  dock  of  the  digital  scope.  This  is  carried  out  to  make  the  number  of 
sampled  phases  of  the  input  signal  to  be  the  smallest.  Once  the  number  of  discrete  phases 
of  the  samples  are  determined,  the  samples  can  then  be  sorted  out  easily  into  groups  of 
identical  input  phases.  For  example,  if  ^he  LO  and  input  signal  are  set  at  1  GHz  and  1020 
MHz  respectively,  the  resultant  baseband  signals  are  at  20  MHz.  When  the  baseband 
signals  are  phas^ocked  to  the  sampling  scope  at  100  MHz,  for  every  cyde  of  the  input 
signal,  exactly  five  samples  are  generated  with  a  phase  difference  of  72  degrees.  These  same 
five  input  phases  of  the  input  signal  are  sampled  repeatedly  for  subsequent  cydes  to 
produce  five  identical  groups  of  samples.  In  the  measurement,  the  phase  of  the  input  signal 
is  kept  phase-4ocked  to  th^  sampling  scope  to  better  than  0.036  degrees  over  a  total 
number  of  5,000  samples.  The  in-i)hase  and  quadrature  signals  are  first  compensated  to 
remove  the  imbalance  errors  and  DC  offsets  by  using  the  procedure  as  outlined  in  Section 
4.0.  A  total  of  5,000  samples,  with  five  groups  of  1,000  samples  each,  are  then  used  to 
compute  the  means  and  RMS  values  of  the  envelope  and  phases  of  the  input  signal.  A  large 
signd  is  first  applied  to  the  scope  to  obtain  approximately  90%  total  deflection.  At  this 
high  SNR,  the  output  amplitude  of  the  envdope  is  approximately  of  Gaussian  distribution 
as  given  by  Eq.(39y  When  the  number  of  samples  is  large,  the  ratio  of  the  measured  RMS 
value  to  the  mean  is  approximatdy  equal  to  (r/a  and  thus  the  input  SNR  is  approximately 

given  by  [a7(2<r’)l.  Next,  the  input  signal  is  attenuated  by  3  dB  and  the  means  and  RMS 
values  of  the  envelope  and  phases  are  then  computed  again. 

Using  the  above  set-up  and  measurement  procedures,  the  statistics  on  the  measured 
envelope  and  phases  are  generated  at  a  baseband  frequency  of  20  MHz  and  40  MHz. 
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As  expected,  the  statistics  are  found  to  be  independent  of  frequency. 

The  RMS  phase  error  in  degrees  as  a  function  of  input  SNR  in  dB  is  plotted  in 
Fig.  12  for  the  case  when  the  baseband  frequency  is  at  40  MHz.  The  circles  are  the 
experimental  points  while  the  solid  curve  is  obtained  by  using  Eqs.(26),  (31)  and  (32).  As 
can  be  seen  from  Fig.  12,  the  experimental  measurement  agrees  very  well  with  theory.  For 
small  input  SNR,  the  RMS  phase  error  is  quite  large.  As  the  input  SNR  is  increased,  the 
RMS  phase  error  decreases  rapidly,  and  for  large  input  SNR,  the  RMS  v^ue  decreases 
inversely  proportional  to  the  square  root  of  the  input  SNR. 

The  ratio  of  the  RMS  value  to  the  mean  of  the  envelope  in  dB,  at  a  baseband 
frequency  of  40  MHz,  versus  input  SNR  in  dB  is  plotted  in  Fig.  13.  The  theoretical  curve  is 
calculated  by  using  Eqs.(33)  ,(35)  and  (36).  As  expected,  the  experimental  measurement 
agrees  very  well  with  the  theory.  As  discussed  in  Section  5.2  ,  the  probability  density 
function  is  a  Rayleigh  distribution  for  small  SNR.  As  a  result,  the  ratio  of  the  standard 

deviation  to  the  mean  is  a  constant  and  equals  to  [(4  —  t)/t  or  —  5.63  dB.  For  large 
input  SNR,  the  ratio  of  the  standard  deviation  to  tne  mean  in  dB  is  inversely  proportional 
to  input  SNR  in  dB. 

Once  the  samples  on  the  phase  of  the  signal  are  generated,  its  approximate 
instantaneous  frequency  can  also  be  determined  by  using  the  relationsMp  given  by  Eq.(7). 
An  experiment  is  carried  out  to  evaluate  the  instantaneous  frequency  error  as  a  function  of 
input  SNR  and  baseband  frequency.  A  band—pass  filter  centered  at  750  MHz  with  a  3-dB 
bandwidth  of  40  MHz  is  inserted  between  the  two  IF  amplifiers  and  the  I/Q  network.  The 
same  procedure  as  described  above  is  used  to  compute  the  phase  values  or  the  input  signal 
as  the  input  SNR  is  varied.  The  LO  frequency  is  set  at  750  MHz  with  an  output  power  of 
10  dBm.  Three  different  IF  frequencies  of  755  MHz,  760  MHz  and  770  MHz  are  used  to 
gaierate  baseband  frequencies  corresponding  to  5  MHz,  10  MHz  and  20  MHz  respectively. 
For  instance,  when  the  frequency  of  the  baseband  signal  is  at  20  MHz,  five  groups  of 
identical  input  phase  samples  of  1000  each  are  generated.  The  RMS  value  on  the  difference 
in  phase  between  consecutive  samples  are  calculated  and  then  normalized  by  the  RMS 
phase  error  of  the  samples.  The  normalized  RMS  phase  difference  as  a  function  of  input 
SNR  is  plotted  in  Fig.  14  for  the  three  different  baseband  frequencies.  The  cross-covariance 
between  consecutive  samples  is  also  computed  and  then  normalized  by  the  RMS  phase 
error  of  all  the  samples.  The  normalized  cross-covariance  as  a  function  of  input  SNR  is 
plotted  in  Fig.  15.  By  comparing  the  two  plots,  there  is  a  close  correspondence  between  the 
two.  In  other  words,  when  the  cross-covariance  is  small,  the  phase  error  from  sample  to 
sample  is  less  correlated  and  thus  a  higher  RMS  phase  difference. 

A  simulation  has  also  been  carried  out  to  verify  the  results  in  Figs.  14  and  15. 

White  Gaussian  random  noise  with  zero  mean  is  first  filtered  by  a  Butterworth  low— pass 
filter  of  20— MHz  bandwidth  before  it  is  added  to  the  baseband  signals.  The  Butterworth 
filter  is  of  1st  order  having  less  than  3  dB  of  attenuation  at  20  MHz  and  15  dB  down  by  30 
MHz.  The  in-phase  and  quadrature  noises  are  uncorrelated.  The  baseband  frequencies  used 
are  DC,  5  MHz,  10  MHz  and  20  MHz.  The  RMS  value  on  the  difference  in  phase  between 
consecutive  samples  is  calculated  and  then  normalized  by  the  RMS  phase  error  of  the 
samples.  The  normalized  RMS  phase  difference  as  a  function  of  input  SNR  is  plotted  in 
Fig.  16. 
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Fig.  13  Ratio  of  RMS  Amplitude  Error  to  Mean  Versus  Input  SNR 


35 


Fig.  14  Normalized  RMS  Phase  Difference  Versus  Input  SNR 

Using  Experimental  Data  (  Noise  Bandwidth  =  20  MHz  ) 
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Fig.  15  Normalized  Cross-Covariance  Versus  Input  SNR 

Using  Experimental  Data  (  Noise  Bandwidth  =  20  MHz  ) 
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Fig.  16  Normalized  RMS  Phase  Difference  Versus  Input  SNR 
Using  Simulated  Data  (  Noise  Bandwidth  =  20  MHz  ) 


At  large  input  SNR,  the  phase  error  due  to  noise  before  it  is  applied  to  the  I/Q 
network  can  m  approximated  by  [10] 

A^^(t)  =  tan-‘[^n8(t)/[a(t)+nc(t)]j  a  n8(t)/a(t)  (41) 

where  the  input  noise  is  represented  by  its  narrow-band  quadrature  components  nc(t)  and 
n8(t).  Both  ncft)  and  n8(t)  are  of  Gaussian  distribution  with  zero  mean  and  of  equal 
variance.  Furtnermore,  they  are  uncorrelated.  Both  nc(t)  and  n8(t)  have  a  spectral  density 
of  ri  over  the  frequency  range  of  |f|  <  As  discussed  in  Section  5.0,  after  both  the 

noise  and  signal  are  mixed  down  to  baseband,  the  phase  error  is  still  approximately  given 
by  Eq.(41).  The  standard  deviation  of  the  phase  error  is  given  by 

"t  =  »/>(«,)  =  1/(2  (42) 


Vq  is  found  to  be  independent  of  the  baseband  frequency  of  the  signal.  When  the 
0 

baseband  frequency  of  the  GW  signal  is  taken  into  account  to  obtain  the  absolute  phase,  a 
linear  phase  progression  is  introduced  as 

^(t)  =  «  fi{\)  +  n8(t)/a(t)  (43) 

where  ^(t)  is  the  phase  function  of  the  baseband  signal  and  is  directly  proportional  to  the 
baseband  frequency.  When  the  baseband  frequency  is  not  DC,  the  spectral  density  of  ns(t) 
is  shifted  by  the  baseband  frequency.  As  a  result,  the  spectral  density  of  the  noise  which  is 
centered  at  DC  has  been  shifted  to  a  higher  frequency  by  the  baseband  frequency  of  the 
signal.  When  the  spectral  density  of  the  noise  is  shifted,  the  cross-correlation  between 
consecutive  phase  samples  is  affected.  In  the  case  where  the  baseband  frequency  is 
comparable  to  the  noise  bandwidth  and  is  much  less  than  the  sampling  frequency,  a  shift  in 
the  noise  spectral  density  by  the  baseband  frequency  is  expected  to  have  a  greater  impact 
on  the  cross— correlation  of  the  phase  samples.  The  phase  error  from  sample  to  sample  is 
expected  to  be  less  correlated,  thus  giving  a  higher  normalized  RMS  phase  difference. 

The  normalized  cross-covariance  is  also  computed  and  plotted  in  Fig.  17.  As 
expected,  there  is  a  close  correspondence  between  the  RMS  phase  difference  and  its 
cross-covariance. 

By  comparing  the  experimental  measurement  to  the  result  of  the  simulation,  there 
is  a  general  agreement  between  the  two.  The  experimental  results  are  less  accmate  at  high 
input  SNR  where  the  phase  errors  due  to  noise  are  very  small  and  the  phase  noise  from  the 
synthesizers  is  significant.  Another  source  of  error  is  that  the  signal  source  and  local 
oscillator  are  not  truly  phase-locked  with  the  dock  of  the  sampling  scope. 
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Fig.  17  Normalized  Cross-Covariance  Versus  Input  SNR 

Using  Simulated  Data  (  Noise  Bandwidth  =  20  MHz  ) 


Another  simulation  has  carried  out  by  increasing  the  bandwidth  of  the  Butterworth 
low— pass  filter  to  40  MHz.  The  Butterworth  filter  is  of  3rd  order  having  less  than  3  dB  of 
attenuation  at  40  MHz  and  15  dB  down  by  45  MHz.  The  sampling  rate  is  still  100  MHz. 
The  normalized  RMS  phase  difference  as  a  function  of  input  SNR  is  plotted  in  Fig.  18.  In 
this  case,  the  noise  spectral  density  is  almost  uniform  across  the  full  frequency  band,  and 
the  consecutive  phase  error  is  existed  to  be  much  less  correlated  as  shown.  As  the 
baseband  frequency  of  the  signal  is  slightly  increased,  the  spectrsd  density  of  the  noise  will 
be  shifted  up  and  aliasing  wm  become  more  severe.  The  exact  cross-correlation  will  be 
more  complicated  to  analyze.  The  normalized  cross-covariance  is  also  computed  and 
plotted  in  Fig.  19. 

Once  the  phase  errors  between  consecutive  samples  are  known,  the  instantaneous 
frequency  error  can  then  be  computed  directly,  when  the  phase  error  from  sample  to 
sample  is  uncorrelated,  the  normalized  RMS  phase  difference  is  expected  to  be  equal  to 

(2)^' In  most  practical  cases,  this  would  represent  the  upper  limit  on  the  error. 


7.3  I/Q  Demodulation  of  Signals 

In  this  section,  two  examples  are  used  to  illustrate  the  effectiveness  of  calibration 
and  moving  averages  on  the  demodulation  of  radar  signals.  Implementation  of  the  I/Q 
network  is  carried  out  by  the  use  of  a  quadrature  IF  mixer  from  Anaren  Inc..  The  LO 
frequency  is  set  at  2  GHz  with  an  output  power  of  10  dBm.  The  in-phase  and  quadrature 
signals  are  digitized  using  the  LeCroy  9450  dual— channel  8— bit  digital  scope.  The  radar 
signals  are  generated  using  the  HP  8791  Frequency  Agile  Signal  Simulator.  The  process  of 
compensation  and  computation  is  carried  out  in  software. 

7.3.1  Pulse  Modulated  CW  Signal 

A  pulse  modulated  CW  signal  at  1960  MHz  with  a  pulse  width  of  approximately  0.5 
^s  is  applied  to  the  I/Q  network.  Figures  20  and  21  show  the  measured  envdope  and 
instantaneous  frequency  respectively  of  the  pulse  modulated  CW  signal.  The  raw  data  of 
the  envelope  is  plotted  in  Fig.20(a)  for  an  input  SNR  of  approximately  28  dB.  The 
imbalance  errors  and  DC  offsets  at  this  frequencv  are  measured  to  be  as  follows;  R  =  0.847, 
(^q  “  ^i)  =  2.92  degrees  ,  aqo/[aKq/2]  =  —  0.55%  and  aip/[aKi/2]  =  5.1  %.  The  DC  offsets 
are  very  small  as  compared  to  the  phase  and  amplitude  imbalances.  As  discussed  in  Section 
3.1,  the  computed  ripples  are  expected  to  have  a  dominant  frequency  of  80  MHz,  which  is 
twice  that  of  the  baseband  frequency  (40  MHz).  However,  since  the  computed  output  is 
sampled  plotted)  at  an  interval  of  10  ns  per  sample  which  corresponds  to  a  sampling  rate 
of  100  MHz,  there  will  be  aliasing  effects  introduced.  As  a  result,  the  adiased  frequency  of 
the  ripples  is  now  at  20  MHz  (100—80)  as  shown.  Figure  20(b)  shows  the  effect  of  applying 
a  7— point  moving  average.  As  expected  the  low— pass  filtering  has  the  effect  of  reducing  the 
amplitude  of  the  ripples  and  at  the  same  time  has  aJso  rounded  off  the  pulse  edges  slightly. 
A  more  logical  approach  is  to  use  the  compensation  technique  first  to  minimize  the  errors 
caused  by  the  imbalance  errors  and  DC  o&ets  and  then  use  a  moving  average,  as  shown  in 
Fig.  20(c).  The  mean  of  the  envelope  is  measured  to  be  —  13.1  dBm  while  the  RMS  value  is 
—  32.56  dBm.  The  measured  peak-to-peak  RMS  value  {20 -logic  [(mean  +  RMS)/(mean  - 
RMS)]}  is  0.098  dB.  The  ratio  of  the  mean  to  the  RMS  value  is  higher  than  the  input  SNR 
because  not  all  of  the  systematic  errors  have  been  removed.  When  a  moving  average  is 
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Fig.  18  Normalized  RMS  Phase  Difference  Versus  Input  SNR 
Using  Simulated  Data  (  Noise  Bandwidth  =  40  MHz  ) 


Pig.  19  Normalized  Cross-Covariance  Versus  Input  SNR 

Using  Simulated  Data  (  Noise  Bandwidth  =  40  MHz  ) 
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(a)  Raw  Data 

Fig.  20  Envelope  of  a  0.5-/iS  Pulse 

(Baseband  Frequency  =  40  MHz  and  T 
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(b)  After  a  7— point  Moving  Average 


Fig.  20  Instantaneous  Frequency  of  a  0.5— /iS  Pulse 

(Baseband  Frequency  =  40  MHz  and  T  =  10 
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(c)  After  Compensation  and  3— point  Moving  Average 

Fig.  20  Envelope  of  a  0.5— /is  Pulse 

(Baseband  Frequency  =  40  MHz  and  T  =  10  ns) 
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(a)  Raw  Data 

Fig.  21  Instantaneous  Frequency  of  a  0.5— /iS  Pulse 
(Baseband  Frequency  =  40  MHz  and  T  = 
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(b)  After  Compensation 

Fig.  21  Instantaneous  Frequency  of  a  0.5— /xs  Pulse 

(Baseband  Frequency  =  40  MHz  and  T  =  10  ns) 
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(c)  After  Compensation  and  3— point  Moving  Average 

Fig.  21  Instantaneous  Frequency  of  a  0.5— ns  Pulse 
(Baseband  Frequency  =  40  MHz  and  T  =1C 


applied  to  the  noise,  the  resultant  noise  RMS  value  is  found  to  decrease  as  expected.  When 
the  moving  average  is  increased  from  one  sample  to  five  contiguous  samples,  the  RMS 

value  of  the  noise  is  found  to  decrease  by  a  factor  of  «  (5)^/^.  On  the  other  hand,  the  signal 
power  is  relatively  unaffected.  As  a  result,  the  output  SNR  improves  as  the  number  of  • 

samples  used  in  the  moving  average  is  increased.  This  happens  because  when  a  moving 
average  is  applied,  the  effective  video  bandwidth  of  the  I/Q  network  is  reduced  by 
approximately  the  square  root  of  the  number  of  samples  used  in  the  moving  average.  As 
long  as  the  video  bandwidth  is  larger  than  the  video  bandwidth  of  the  signd,  the  measured 
envelope  of  the  signal  is  not  affected. 

Figure  21  shows  the  instantaneous  frequency  distribution.  The  raw  data  shows  that 
a  peak— to— peak  frequency  deviation  is  about  4  MHz.  When  a  3— point  moving  average  is 
used  after  compensation,  the  RMS  instantaneous  frequency  error  is  reduced  drastically  to 
below  80  KHz.  Other  types  of  low— pass  digital  filters  such  as  Butterworth  and  Chebyshev 
are  also  used  on  the  compensated  instantaneous  frequency  of  the  signal.  However,  they  are 
found  to  introduce  ringing  and  overshoot  which  significantly  distort  the  instantaneous 
frequency. 

7.3.2  Linear  FM  Signal 

A  linear  FM  signal  with  a  center  frequency  of  1960  MHz,  and  a  frequency  excursion 
of  10  MHz  is  applied  to  the  I/Q  network.  The  envelope  and  instantaneous  frequency  are 
shown  in  Figs.  22  and  23  respectively.  It  has  the  same  pulse  width  and  input  SNR  as  for 
the  pulsed  signal.  Even  though  the  FM  signal  has  a  frequency  excursion  of  10  MHz,  the 
same  set  of  compensation  parameters  are  used  as  if  it  were  a  CW  signal  centered  at  1960 
MHz.  As  can  be  seen  from  the  plots,  the  compensated  envelope  and  instantaneous 
frequency  of  the  input  signal  have  b^n  reproduced  accurately.  In  other  words,  the 
imbalance  errors  and  DC  of&ets  of  the  I/Q  network  are  relatively  insensitive  to  frequency 
variations.  As  a  result,  coarser  frequency  calibration  steps  can  be  used. 

8.0  SUMMARY  AND  CONCLUSIONS 

A  simple  in— phase/quadrature  (I/Q)  demodulator  architecture  which  can  measure 
the  amplitude,  phase  and  instantaneous  frequency  of  radar  signals  has  been  examined  in 
this  report.  By  splitting  the  in-phase  and  quadrature  components  of  the  input  signal  using 
analog  components,  this  architecture  is  potentially  useful  for  radar  ESM  applications  where 
wide  instantaneous  bandwidth  and  simple  algorithms  for  extracting  the  modulation 
characteristics  of  radar  signals  are  required.  However,  there  are  amplitude  and  phase 
imbalance  errors  and  DC  offsets  in  practical  I/Q  networks,  which  can  introduce  large 
systematic  errors  to  the  measurement.  Two  techniques,  compensation  and  low— pass 
filtering  using  a  moving  average,  have  been  proposed  in  this  report  for  improving  this 
simple  architecture. 

The  effect  of  amplitude  and  phase  imbalances  and  DC  offsets  is  found  to  produce 
ripples  on  the  measurement;  the  frequency  of  the  ripples  is  harmonically  related  to  the 
baseband  frequency  of  the  signal.  For  a  given  set  of  imbalance  errors  and  DC  offsets,  the 
error  introduced  by  the  ripples  on  the  instantaneous  frequency  is  more  severe  and  is 
directly  proportional  to  the  baseband  frequency  while  the  error  is  constant  for  both  the 
amplitude  and  phase. 
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Fig,  22  Envelope  of  a  Linear  FM  Pulse 

(  Pulse  Width  =  0.5  (is,  Center  Frequency  =  40  MHz, 
Af  =  10  MHz  and  T„  =  10  ns) 
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(b)  After  Compensation  and  5— point  Moving  Average 


Pig.  22  Envelope  of  a  Linear  FM  Pulse 

(  Pulse  Width  =  0.5  ^s,  Center  Frequency  =  40  MHz, 
nf  =  10  MHz  and  T  =  10  ns) 
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(b)  After  CompensatioTi  and  5-point  Moving  Average 


Fig.  23  Instantaneous  Frequency  of  a  Linear  FM  Pulse 

(  Pulse  Width  =  0.5  fis,  Center  Frequency  =  40  MHz, 
af  =  10  MHz  and  T  =  10  ns) 
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Some  I/Q  networks  have  been  evaluated  and  the  mean  imbalance  errors  and  DC 
onsets  over  a  large  frequency  range  can  be  large.  In  order  to  keep  the  systematic  errors 
small,  some  form  of  compensation  is  needed  to  remove  the  mean  errors  and  offsets.  Once 
the  mean  errors  and  of&ets  are  eliminated,  the  residual  RMS  errors  as  a  function  of 
frequency  and  input  power  level  are  quite  small.  If  further  improvement  is  required,  a 
calibration  dependent  on  frequency  and  input  power  level  may  be  needed. 

The  measurement  accuracy  on  the  amplitude,  phase  and  instantaneous  frequency 
has  also  been  analyzed  and  given  in  terms  of  input  SNR.  For  large  inpu*  SNR,  the  standard 
deviation  of  the  phase  error  decreases  inversely  proportional  to  the  square  root  of  the  input 
SNR  while  the  standard  deviation  of  the  envelope  remains  constant  as  the  input  signal 
level  is  varied. 

There  are  two  basic  functions  performed  by  using  a  moving  average.  The  effective 
noise  bandwidth  of  the  I/Q  network  can  be  reduced  with  no  appreciable  degradations  to 
the  signal  of  interest,  with  consequent  improvement  of  the  output  SNR.  The  ripples 
produced  by  the  imbalance  errors  and  DC  offsets  are  usually  outside  the  video  bandwidth 
of  the  signal  of  interest  and  thus  can  be  reduced  by  low— pass  filtering.  In  addition,  a 
moving  average  is  easy  to  implement  digitally  and  does  not  introduce  distortions  such  as 
ringing  and  overshoot. 

Both  the  techniques  of  compensation  and  low— pass  filtering  using  a  moving  average 
have  been  successfully  demonstrated  on  the  demodulation  of  pulsed  and  linear  FM  signals. 
They  have  been  shown  to  be  effective  for  improving  the  output  signal— to— noise  ratio  and 
reducing  systematic  errors. 
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